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Abstract 

Local heat transfer coefficients were measured 
on a smooth and roughened NACA 0012 airfoil. 
Heat transfer measurements on the 0.533 m chord 
airfoil were made both in flight on the NASA Lewis 
Twin Otter Icing Research Aircraft and in the 
NASA Lewis Icing Research Tunnel (IRT). 
Roughness was obtained by the attachment of uni- 
form 2 mm diameter hemispheres to the airfoil 
surface in four distinct patterns. Flight data were 
taken for the smooth and roughened airfoil at 
various Reynolds numbers based on chord in the 
range 1.24xl0 6 to 2,50xl0 6 and at various angles of 
attack up to 4°. During these flight tests, the free 
stream velocity turbulence intensity was found to 
be very low (<0.1%). Wind tunnel data were ac- 
quired in the Reynolds number range 1.20X10 6 to 
4.52xl0 6 and at angles of attack from -4° to +8°. 
The turbulence intensity in the IRT was 0.5 to 0.7% 
with the cloud-generating sprays off. 

A direct comparison was made between the re- 
sults obtained in flight and in the Icing Research 
Tunnel. The higher level of turbulence in the IRT 
versus flight had little effect on the heat transfer for 
the lower Reynolds numbers but caused a moderate 
increase in heat transfer at the higher Reynolds 
numbers. Turning on the cloud-generating spray 
nozzle atomizing air in the IRT did not alter the 
heat transfer. Roughness generally increased the 
heat transfer by locally disturbing the boundary 
layer flow. Finally, the present data were compared 
with leading edge cylinder and flat plate heat 
transfer values which are often used to estimate 
airfoil heat transfer in computer codes. 


Nomenclature 

A surface area of gage 
c chord length 

Fr Frossling number based on chord 

h local heat transfer coefficient 

k thermal conductivity of air 

M Mach number 

Nu Nusselt number based on chord 

Q et electric power input to heater 

Q eTtd heat loss from unguarded end of heater 

Q^p heat loss through gap 

Q rad heat loss due to radiation 

Re Reynolds number based on chord 

T f film temperature 

T # measured static temperature 

T t total temperature 

T g measured gage temperature 

V free stream velocity 

Greek symbols 

e surface emissivity of polished aluminum 
p viscosity of air 
p density of air 
a Stefan-Boltzman constant 

Introduction 

The problem of ice accretion on aircraft has 
historically drawn considerable attention. The 
hazards associated with ice formation on wings and 
engine inlets are well studied and quite apparent. 
Glaze or rime ice formations on an airfoil add 
weight, reduce lift, increase drag and may cause 
flight control problems [1-4]. A further danger 


arises when pieces of ice shed off a wing or engine 
inlet and cause damage to parts of the aircraft such 
as the fan and compressor blades, tail, rotors, etc. 
Ice formation on rotary aircraft, has drawn much 
study because of the difficulty of preventing or 
removing it [5]. These problems have once again 
come into focus for civilian and military aircraft, 
especially in the application of low-flying missiles 
and helicopters in regions of cooler climate. 

The two general methods used to alleviate the 
problems of ice accretion are anti-icing and de-icing. 
Anti-icing methods prevent ice from forming, most 
often by heating the affected area above the water 
freezing temperature. De-icing methods, on the 
other hand, remove ice after it begins to grow but 
before it causes much adverse effect. This is gener- 
ally accomplished by melting or mechanically 
cracking the ice and allowing centrifugal or 
aerodynamic forces to shed it. In order to most 
efficiently apply either of these methods, it is first 
necessary to attain a good understanding of the 
icing phenomenon and to be able to predict whether 
or not ice will grow under specific environmental 
conditions and for specific locations. Further, ice 
growth prediction is especially vital in applications 
in which no ice protection equipment is used. Key 
to this prediction is the heat transfer. 

If an aircraft passes through cool air contain- 
ing supercooled liquid water droplets, an energy 
balance shows that if the convective and 
evaporative cooling, as well as the warming of the 
impinging supercooled droplets to the freezing tem- 
perature, can overcome the kinetic and viscous 
heating and thus remove the latent heat of fusion 
from the impinging droplets, ice will form on the 
surface [6,7,8]. The dominant heat loss terms in the 
thermal analysis are convective and evaporative 
cooling; hence the importance of the local convective 
heat transfer coefficient, h. 

Icing facilities and ice accretion modeling 
codes must accurately simulate heat transfer in 
natural icing conditions. Heat transfer is 
dependent on the free stream turbulence level, and 
since wind tunnels typically have higher turbulence 
levels than are found in flight, there may be a prob- 
lem in using wind tunnels to predict flight condi- 
tions. In this study turbulence intensity has been 
measured up to a free stream velocity of 94 m/sec 
(210 miles per hour) to be equal to or less than 0.7% 
in the NASA Lewis Icing Research Tunnel (IRT). 
Turbulence intensity for flight conditions measured 


during this study have been found to be too low to 
make a meaningful measurement (<0.1%) for 
smooth air. Previous measurements in rough air 
below a layer of cumulus clouds have given 
somewhat higher levels (0.2-0.4%). 

One objective of the present tests was to 
determine the differences in local heat transfer from 
a smooth and roughened airfoil between flight and 
tunnel conditions. A second objective of this work 
was to obtain much needed airfoil heat transfer 
data to better describe the thermal physics 
occurring during the icing process and, specifically, 
to provide accurate airfoil heat transfer data for use 
in ice growth prediction computer codes. One such 
code, NASA’s LEWICE [9], currently incorporates 
an integral boundary layer subroutine to calculate 
heat transfer coefficients. Often heat transfer from 
an airfoil is estimated with cylinder-in-crossflow 
heat transfer data for the stagnation region and flat 
plate heat transfer coefficients for the rest of the 
airfoil surface. The present tests provide actual 
data for a NACA 0012 airfoil for a smooth surface 
as well as for quantifiable roughness patterns. 

Limited data currently exist for heat transfer 
from an airfoil. A NACA study (1946-1951) [10,11] 
compared in-flight heat transfer from an airfoil, in 
clear air and during icing conditions, with results 
from the IRT. For the flight data two separate 
airfoils, a NACA 0012 and a NACA 65,2-016, were 
tested at a 0 degree angle of attack, while only the 
65,2-016 was subsequently used in the IRT. In the 
"flat plate" region (i.e., the region away from the 
stagnation area), the data showed a substantial 
difference between flight and IRT heat transfer on 
the forward portion of the airfoil where the 

boundary layer was laminar. The IRT data was 
over 30% higher than the flight data. This 
difference has been attributed to the higher 

turbulence intensities present in the IRT. This 

conclusion is also supported by the fact that the 
flight and IRT data agreed fairly well on the 
downstream portion of the airfoil where the 

boundary layer was assumed fully turbulent. 

Besides being restricted to a 0° angle of 
attack, two other factors limit the usefulness of the 
previous data for computer code predictions. First, 
the data set is incomplete and somewhat 
inconsistent in the stagnation region, the area 
where ice growth initiates. Secondly, data were not 
taken for a rough surface, which can significantly 
alter boundary layer characteristics and thus the 
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local heat transfer. Roughness, the result of early 
ice growth, may force a laminar boundary layer into 
transition in the ice formation zone. This behavior 
was observed in recent experiments performed on a 
cylinder in crossflow under different turbulence and 
roughness conditions [12]. 

A 1985 study at the University of Kentucky 
[13] included the determination of heat transfer 
coefficients from a smooth NACA 0012 airfoil in a 
subsonic wind tunnel, as well as from a five minute 
ice accretion shape. The smooth airfoil 
measurements were taken at various angles of 
attack (-8° through +8°) and for a chord-based 
Reynolds number range of 7.6x10 s to 2.0xl0 fi . The 
0° angle of attack data generally agreed with the 
NACA study and the data showed a larger angle 
dependence on the suction side as compared to the 
pressure side. Again, however, more complete 
roughness data is lacking. 

The present study focused on heat transfer 
measurements on a NACA 0012 airfoil. The NACA 
0012 was chosen because it is a symmetric profile 
that is commonly used in helicopter main rotor and 
tail rotor applications where it may not be possible 
to control ice growth by electric heating or 
pneumatic boots. Local heat transfer coefficients 
were calculated from measurements taken on the 
forward 8 percent chord of a smooth and roughened 
NACA 0012 airfoil with a 0.533 m (21 in) chord 
length. Roughness was obtained by the attachment 
of small hemispheres of uniform size (2 mm diame- 
ter) onto the airfoil in a set and reproducible 
pattern. Several distinct position patterns, similar 
to those employed by Schlicting [14] in his boundary 
layer work, were used. These patterns were chosen 
to facilitate numerical modeling of the roughness in 
various computer codes. Heat transfer measure- 
ments were recorded in flight on the NASA Lewis 
Twin Otter icing research aircraft and in the NASA 
Lewis Icing Research Tunnel. Flight data were col- 
lected for smooth and roughened surfaces at various 
aircraft speeds and various angles of attack up to 
4°. Data were acquired in the IRT for smooth and 
roughened airfoil surfaces at various tunnel 
airspeeds, with and without spray nozzle atomizing 
air, and for various angles of attack from -4° (heat 
flux gages on the pressure side) to +8° (heat flux 
gages on the suction side). Detailed results of the 
flight tests, including comparisons with published 
data, have previously been reported in Reference 
15. Selected data from the IRT and flight tests will 
be presented here as Frossling number versus 


position on the airfoil for various roughness 
patterns and angles of attack. Stagnation region 
data are compared with Frossling's cylinder-in- 
crossflow solution [16] and data further aft on the 
airfoil are compared with flat plate correlations 
[17]. Complete results are available in Reference 
18. 

Experiment Equipment and Procedure 
Apparatus 

UACA QQ12 Test Airfoil 

Heat transfer measurements were made on a 
NACA 0012 airfoil which was made of mahogany 
and had a chord length of 0.533 m (21 in) and a 
span length of 1.8 m (6 ft). The airfoil is shown in 
Figure 1. 

An array of heat transfer gages was located in 
a removable section at the center of the span. Each 
heat flux gage was equipped with a thermocouple 
and foil heater which were connected to individual 
circuits that allowed the gages to be operated in a 
constant temperature mode. Figure 2 shows a 
cross-section of the airfoil and the location of the 
heat transfer gages. Table 1 lists the surface 
distance from the geometric stagnation point to the 
center of each gage and its heat transfer surface 
area. The airfoil actually contained twenty-eight 
heat flux gages but only twelve of them in the 
forward region were used in this study. A comput- 
erized data acquisition system was employed to 
record voltage and current input to each gage as 
well as gage temperature. Details of the model and 
instrumentation can be found in References 15 and 
18. 

Surface roughness was added to the airfoil by 
fastening hemispheres of silver alloy to the surface 
with cyanoacrylic adhesive. The hemispheres were 
2 mm in diameter and were attached to the airfoil 
in different patterns. Three patterns were 
employed in the flight tests while four were used in 
the tunnel tests. Figure 3 shows planar sketches of 
the location of the roughness elements relative to 
the heat flux gages for each of the patterns. The 
thermal resistance of the gage surface was not 
altered significantly because of the sparse spacing 
of the elements and the high conductivity of the 
silver alloy. A numerical heat conduction computer 
model predicted that the temperature at the tip of 
each roughness element was less than 0.28°C 
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(0.5°F) lower than the temperature of the smooth 
aluminum gage surface. Also, each roughness 
element generally increased the heat transfer 
surface area by only one percent per gage. No 
attempt was made to account for the presence of the 
roughness elements in the data reduction. 

Twin Otter Aircraft 

The NACA 0012 airfoil was flown atop the 
NASA Lewis Twin Otter Icing Research Aircraft. 
Details of the aircraft and parameter 
measurements are found in References 15 and 18. 
With the NACA 0012 research airfoil mounted 
vertically on the top of the aircraft the maximum 
sustainable speed was around 69 m/sec (154 miles 
per hour) at 1585-2250 m (5200-7400 ft) pressure 
altitude and a temperature range of 10-2 1°C (50- 
70°F). Variance of the test airfoil angle of attack 
was achieved through yawing the aircraft with a 
combination of rudder and aileron. Flight tests 
were performed in darkness to avoid any effect that 
solar heating may have had on the heat flux gage 
temperatures. 

Icing R esearch . T unnel 

The NACA 0012 heat transfer airfoil was also 
tested in the NASA-Lewis Icing Research Tunnel 
(IRT). The IRT is a closed-loop subsonic 
refrigerated air tunnel used primarily for icing 
studies. A plan view of the IRT is shown in Figure 
4. Air flow is induced by a 5000 horsepower, 7.31 m 
(24 ft) diameter fan, and airspeed in the 1.83 by 
2.74 m (6 by 9 ft) test section can be varied from 
about 9 to 125 m/sec (20 to 280 miles per hour). The 
air is cooled by passage through a heat exchanger 
unit which maintained total air temperature 
around -6.7°C (20°F), plus or minus 0.55°C (1°F). 
The tunnel air temperature was measured by an 
array of eleven thermocouples positioned on the 
first turning vane upstream of the test section. 
Spray bars, located approximately 14.6 m (48 ft) 
upstream of the test section, are used to produce the 
icing cloud. The spray bars contain a collection of 
spray nozzles which shoot a combination of 
pressurized air and water to yield a continuous and 
uniform cloud of very small supercooled water 
droplets. The spray bars are also heated with a 
separate closed steam loop to prevent nozzle freeze- 
up. The liquid water content of the tunnel icing 
cloud can generally range from about 0.2 to 3.0 
gr/m 3 , and the drop diameter can range from about 
5 to 40 pm [19]. Typical cloud conditions require 


the nozzle atomizing spray air to be set at 60 psi 
and due to facility constraints roughly 82°C (180°F). 
In the present heat transfer tests, ice growth on the 
airfoil was to be avoided; therefore, when the spray 
bars were employed, only the 60 psi air was used. 
No water was passed through the nozzles. 

Airspeed was measured using a Pitot-static 
probe attached to the wall in the tunnel test section. 
It was positioned roughly 3.05 m (10 ft) in front of 
the test airfoil, about 63 cm (25 in) from the ceiling 
and 51 cm (20 in) from the wall. The velocity 
determined at this point was assumed to be the free 
stream velocity at the center of the tunnel and 
based on many years experience with the tunnel 
this assumption is correct within 0.9-1. 3 m/sec (2-3 
miles per hour). 

The airfoil was mounted vertically on the 
turntable in the floor of the test section as shown in 
Figure 5. The lower end was fastened to the floor 
turntable so that the chord of the airfoil was 
parallel to the tunnel walls when the turntable 
angle was set at zero degrees. The top of the airfoil 
was fitted with a metal collar and attached to the 
tunnel ceiling at a single point such that the airfoil 
could rotate with the turntable. 


Turbulence Me asurements 

Turbulence measurements were made with a 
standard constant temperature hot wire 
anemometer operating in an uncalibrated mode 
[20]. Turbulence measurements for the flight tests 
are described in References 15 and 18. 

Turbulence measurements were made in the 
IRT with the same constant temperature hot wire 
equipment that was used in the flight tests. The 
probe, consisting of a single tungsten wire, was 
positioned in the center of the tunnel about a foot in 
front of the gages on the test airfoil. Turbulence in- 
tensities were measured with the cloud making noz- 
zle atomizing air sprays off. Turbulence measure- 
ments with the cloud making air sprays tu rne d on 
were found to be difficult to measure with hot wire 
equipment because the ambient tunnel air 
temperature generally averaged 6.7°C (20°F) while 
the spray air temperature had to be maintained 
near 82°C (180°F) due to facility limitations. Since 
hot wire equipment is a heat transfer sensing 
device, turbulence readings would be affected by 
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this temperature difference and would not measure 
true turbulent velocity fluctuations. 

Heat Transfer Measurements. 
flight Measurements 

All heat transfer data acquisition flights were 
made in darkness to avoid solar radiation on the 
gages and the airfoil. Flights were conducted at an 
altitude that provided smooth atmospheric 
conditions. 

When steady-state air flow conditions were 
established, the heated aluminum strips were ad- 
justed to a uniform temperature which was 
typically in the range of 32-4 1°C (90-105°F). When 
steady state thermal conditions were reached, data 
recording was initiated. About two minutes were 
required to obtain and record the required ten scans 
of all data channels. 

Icing Research Tunnel Measu re ments. 

The heat transfer tests done in the IRT were 
performed in much the same way as the flight tests. 
After the initial tunnel cooldown, the heat transfer 
experiments were begun. The airfoil angle of attack 
was set by rotating the turntable to the desired 
position and the tunnel air velocity was varied by 
adjusting the fan speed. When steady state tunnel 
conditions were achieved, the airfoil heaters were 
adjusted to a uniform temperature, typically in the 
range of 32 to 38°C (90 to 100°F). Again, data 
recording was initiated after steady state thermal 
conditions were obtained. Roughly two minutes 
were required to obtain and record ten scans of all 
data channels. Runs were made with and without 
the cloud generating spray air turned on. 

Data Reduction 

The average heat transfer coefficient from 
each gage was obtained from the applied voltage 
and current and the calculated temperature 
difference between the gage and the free stream 
total temperature. Since only the convective heat 
transfer was desired, the radiation heat loss had to 
be subtracted from the total electric power input to 
each heater. Further, the heater gages embedded 
in the airfoil were secured in place and separated 
from each other by an epoxy resin. Some heat was 
conducted from the edges of each heater gage 
through the epoxy and convected from the surface of 


the airfoil in the gaps between the gages and from 
the unguarded ends of the gages. These losses were 
also subtracted from the electric power. Therefore, 
the local heat transfer coefficient, h, for each 
aluminum heater gage was calculated from: 

h=(Q E ,-Qr*i-Qg Bp -Qend) ' (ACT/T,)) (1) 

where Q EI (voltage x current) is the total electric 
power input to each heater. The quantity is the 
radiation heat loss, which is estimated by: 

Q rad = a A e (T r 4 - T t 4 ) (2) 

A value of 0.045 was used for e, the emissivity of 
polished aluminum, and a is the Stefan-Boltzman 
constant. The quantities Q^ p and Q end are the heat 
lost through the epoxy gaps separating the 
aluminum gages and the unguarded ends of the 
heaters, respectively. These quantities were 
obtained from an exact solution for heat conduction 
in a rectangle with appropriate boundary 
conditions, which is detailed in Reference 18. The 
remaining quantities are: A, the surface area of 
each aluminum strip; T , the measured wall 
temperature of each gage; and T t , the total 
temperature. For the flight data, T t was calculated 
from the measured static temperature, T a , and the 
true airspeed, i.e., 

T t = T a (l + M/5) (3) 

where M is the Mach number. Typical values for T t 
ranged from 10 to 21°C (50 to 70°F). For the IRT 
data, T t was measured with two thermocouples 
positioned on the leading edge of the airfoil and was 
typically around -6.7°C (20°F). 

The heat transfer data in this analysis is pre- 
sented as Frflssling number based on chord length 
which is calculated as: 

Fr = Nu / Re 06 = [(hc/k) / (pVc/^i) 05 ] (4) 

where c is the 0.533 m (21 in) chord length. The 
density, p, was calculated from the ideal gas 
relation for air using the static temperature and 
pressure at the test airfoil location. The velocity, V, 
for flight data was measured with Twin Otter 
instrumentation and converted to the true test 
section velocity with the aid of calibrations obtained 
on previous Twin Otter flights, while the velocity 
for the tunnel data was measured with the IRT 
Pitot static probe. The thermal conductivity, k, and 
viscosity, p, were obtained as functions of 
temperature from air data in Reference 21. These 
thermal properties were evaluated at the film 
temperature given by: 

T f = (T g + T t )/2 (5) 


5 


An uncertainty analysis according to the 
method of Kline and MeClintock [22] was performed 
on the calculated local heat transfer coefficient and 
the Frossling number. The uncertainties for each 
gage were similar and averaged around 4.5% for h 
and 5% for Fr. A substantial part of this uncertain- 
ty was found to be due to uncertainty in the gap 
heat loss term because the thermal conductivity of 
the epoxy gap material was not known exactly and 
was assumed to be 0.11 (with a ±50% uncertainty), 
a typical value for epoxy of this nature. This would 
not be a random error but would tend to bias the 
data either high or low. However, good general 
agreement with flat plate data seemed to confirm 
that the epoxy thermal conductivity value used was 
correct. 

Besets and Discussion 

In this section selected results of the heat 
transfer tests conducted in the Icing Research 
Tunnel (IRT) as well as data obtained during the 
Twin Otter aircraft flights will be presented. 
Special emphasis is given to the effect on heat 
transfer of the different free stream turbulence 
levels. The turbulence intensity measured with hot 
wire equipment during the flight runs was found to 
be veiy low (<0.1%). However, in the IRT with the 
cloud making sprays off, the turbulence intensity 
level was found to be 0.6, 0.52 and 0.7% at tunnel 
air speeds of 31, 63, and 94 m/sec (70, 140 and 210 
miles per hour), respectively. Previous studies 
measured the IRT turbulence levels to be around 
0.5% with the cloud making air sprays off and 
around 2% with spray equipment operating; howev- 
er, this latter result is somewhat suspect because of 
the aforementioned concerns regarding the required 
spray air temperature. 

Generally the IRT data exhibited the same 
trends as the flight data previously reported in 
Reference 15. Figure 6 shows Frossling number 
based on chord as a function of s/c for the smooth 
airfoil at a 0° angle of attack without and with tun- 
nel spray air for several Reynolds numbers. The 
data plotted in this manner collapse onto a single 
curve which shows that the Nusselt number is 
proportional to the square root of the Reynolds 
number. A least squares curve fit of the equation 

Nu=A(ReP (6) 

for the tunnel data [18] showed that the Nusselt 
number did correlate with Reynolds number raised 
to the 0.5 power. It should be noted here that the 


flight Nusselt number data more accurately 
correlated with Reynolds number raised to the 
power of around 0.43. However a sufficiently good 
correlation with the square root of the Reynolds 
number was found to justify presentation of the 
data as Frossling number. The solid line on the fig- 
ure represents the averaged, smooth-surface, 0° 
angle of attack, flight data and will be reproduced 
on subsequent figures for reference. As expected, 
the Frossling number is greatest at the stagnation 
point, with an average value of 4.56, and then trails 
off smoothly to an average value of 1.06 at an s/c 
value of 0.083, The IRT stagnation point average 
value of 4.56 is only 6% higher than the average 
stagnation value of 4.30 for the flight data. The cor- 
relations developed by Lowery and Vachon [23] pre- 
dict a stagnation point heat transfer enhancement 
due to a 0.7% free stream turbulence of 3.5% at 
Reynolds number of l,2x!0 6 and 5.9% at Reynolds 
number of 3.5xl0 6 . The reason for the ,r bump" at s/c 
= 0.048 is unknown; there are no obvious rough 
spots or steps on the surface of the airfoil at this or 
any other point. This same bump was present in 
the flight data. 

It is possible to define the Frossling number 
in terms of the diameter of a cylinder inscribed in 
the leading edge of the airfoil. This allows 
comparison of Frossling's analytical solution for 
heat transfer in the stagnation region of a circular 
cylinder [16] with the present data. For a NACA 
0012 airfoil, the leading edge equivalent diameter is 
3.16% of the chord or 0.0169m (0.664 in) for the 
airfoil tested [24]. The experimental average 
Frossling number based on leading edge equivalent 
diameter for the smooth airfoil was found to be 
0.806, roughly 15% lower than the 0.945 value 
predicted by Frosslings cylinder solution. 
Frossling's analytical results are often used with an 
equivalent leading edge diameter to compute heat 
transfer in the stagnation region for airfoils and 
turbine blades. However, it is uncertain as to 
whether the validity of this method has been proven 
experimentally. 

For comparison purposes, representative 
cases of smooth and roughened tunnel and flight 
data under similar conditions (i.e. same angle of 
attack and comparable Reynolds numbers) will be 
presented. 

Figure 7 shows the heat transfer results for 
both spray conditions of the IRT tests as well as for 
the flight test. This data is for the smooth airfoil at 
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a 0° angle of attack and a nominal Reynolds num- 
ber equal to 1. 2x10 s . Compared with the flight 
data, gages 2-6 show a Frossling number only 2-3% 
higher for the tunnel cases without spray air, and 3- 
5% higher with spray air. Further aft on the airfoil 
(s/c > 0.03), the Frossling number averages 6% 
higher both with and without spray air cases. All of 
these values are within the limits of the calculated 
experimental uncertainty. 

Figure 8 shows similar Frossling number 
behavior for a nominal Reynolds number of 2. 4x10 s . 
The two tunnel cases, with and without spray air, 
again agree quite well (within 2%). However, 
comparing both cases with the flight data shows 
that the Frossling number tunnel data is roughly 
7% higher on gages 2-6 and 10% higher on gages 
greater than 6. 

The 2° and 4° angle of attack smooth-airfoil 
data exhibited similar behavior, showing good 
general agreement between tunnel and flight data 
at the lower Reynolds numbers, while at the higher 
Reynolds numbers the Frossling numbers for the 
tunnel cases were somewhat higher than the flight 
data. 

It seems evident that the small increase in 
turbulence level has a slight affect on the heat 
transfer and this effect is magnified with increasing 
Reynolds number. A consequence of this behavior is 
seen in the Nusselt versus Reynolds number power 
law constants. Generally the exponents on the 
downstream gages of the IRT data are slightly 
higher than the corresponding flight data. Further, 
and perhaps more importantly, the addition of 
spray air (recall that in these experiments only 
spray air and not any spray water was turned on) to 
the tunnel stream does not affect the heat transfer 
and thus apparently does not affect the turbulence 
level. It should be mentioned, however, that it is 
possible that the spray air does increase the tunnel 
turbulence but the leading edge heat transfer is not 
sensitive to this change. 

The addition of surface roughness to the 
leading edge of the airfoil as shown in Figure 3a, for 
0° angle of attack, locally increased the stagnation 
point heat transfer about 10% but did not appear to 
affect the downstream heat transfer. This increase 
can be partially explained by the 4% increase in 
surface area due to the presence of the 
hemispherical roughness elements. It could also be 
attributed to a disturbance of the boundary layer by 


the relatively huge roughness elements followed by 
a return to laminar flow sufficiently past the 
leading edge trip point. 

The sparsely roughened airfoil (Figure 3b) ex- 
hibited similar behavior. The flight and tunnel 
data for the sparsely roughened airfoil at a 0° angle 
of attack and Reynolds numbers 1.2x10 s and 
2.4x10 s are presented in Figures 9 and 10, respec- 
tively. The axial positions of the roughness ele- 
ments are denoted by the arrows above the abscis- 
sa. Compared with the smooth-airfoil cases the 
data show an increase in heat transfer on the 
stagnation gage which is similar to the leading edge 
roughness case. Downstream of the stagnation 
region, the Frossling number exhibits a pattern of 
increasing at and immediately downstream of the 
roughness row position, then falls off slightly. Note 
that the downstream gages show somewhat higher 
Frossling numbers for the higher Reynolds number 
(Figure 10) versus the lower Reynolds number case 
(Figure 11). The power law curve fits showed that 
at an s/c greater than 0.02, the Nusselt number 
begins moving to a Reynolds number dependence 
higher than Re 05 . This trend was especially true for 
the tunnel data and it is an indication of boundary 
layer transition. With few exceptions the flight and 
tunnel Frossling number values on most gages 
agreed within the uncertainty limits for the low 
Reynolds number case. At the higher Reynolds 
number, as with the smooth case, the two tunnel 
conditions agreed quite well, but most gages for 
both cases show Frossling number values 4-10% 
higher than corresponding flight values. Similar 
trends were also observed on the 2° and 4° angle of 
attack data, though the 4° case exhibited somewhat 
higher magnitude; some gages for the 2.4xl0 6 
Reynolds number show tunnel heat transfer data 
up to 16% higher than the flight data. 

Increasing the density of the roughness 
elements from the sparse to the dense 1 pattern 
(Figure 3c) has a dramatic effect on the heat 
transfer downstream. The data for the 0° angle of 
attack, dense 1 roughened airfoil are illustrated in 
Figures 11 and 12. The heat transfer from gages 2- 
8 is substantially higher than the corresponding 
sparsely roughened case. Further downstream past 
gage 7, the density of roughness elements 
decreases, and at gages 10 and 11, the effect of the 
increased density of the roughness elements 
upstream seems to have nearly damped out. This 
trend indicates that if there is roughness of 
sufficient magnitude present, the boundary layer is 
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perturbed locally and immediately downstream. 
However, as the density of roughness is reduced in 
the downstream direction, the heat transfer 
recovers to a level that is consistent with the sparse 
roughness pattern. 

For the lower Reynolds number of 1.2x10*, 
both tunnel data sets again agree; however, they 
show values 2-8% higher than the flight data. The 
2,4x10* Reynolds number case illustrates similar 
behavior, with the tunnel data being 18-25% higher 
than the flight data. This result is somewhat 
surprising because it would seem that as the 
boundary layer is more and more disturbed by the 
roughness, the free stream turbulence would have 
less and less effect on the Frossling number. 
However, in this dense 1 case the free stream 
turbulence appears to have a greater effect on the 
heat transfer than in the sparse roughness cases. 
Moving to a higher angle of attack did not alter this 
trend. 

Frossling data versus s/c for the 0°, dense 2 
roughness pattern (Figure 3d) are presented for low 
and high Reynolds number cases in Figures 13 and 
14, respectively. Note that these graphs show only 
tunnel data, as no dense 2 flight data were ac- 
quired. Examination of Figure 3 shows that the 
dense 2 roughness pattern is very similar to the 
dense 1 roughness pattern, except that the 
roughness element density is constant throughout 
and does not decrease at s/c greater than 0.04 as the 
dense 1 pattern does. As expected, the heat 
transfer results shown in these dense 2 graphs are 
very similar to the results of the dense 1 cases. The 
only notable difference between the dense 2 and the 
dense 1 cases occurs in the region for s/c greater 
than 0.04. The Frossling number in the dense 2 
cases in this region, while still falling off from the 
gage 7 maximum and then leveling off, is somewhat 
higher than the dense 1 data. Thus, the dense 2 
roughness data exhibit the same heat transfer 
behavior as the dense 1 cases, except that the dense 
2 data yield a higher Frossling number in the 
region of higher roughness density (s/c >0.04). 
There is no upstream heat transfer effect due to the 
increase of downstream roughness density. 

Generally, it would seem that the IRT is a 
relatively clean wind tunnel, at least with respect to 
heat transfer in the forward portion of an airfoil. 
The slightly higher turbulence level in the tunnel 
had minimal effect on the heat transfer from the 
smooth airfoil at the lower Reynolds number; the 


higher Reynolds number did illustrate a greater 
turbulence induced heat transfer increase. The 
addition of surface roughness, however, seemed to 
magnify the effect of turbulence on heat transfer. 
The addition of spray atomizing air had virtually no 
effect on the heat transfer, and it is therefore 
believed that the spray air did not significantly 
alter the free stream turbulence level in the tunnel. 

The general Frossling number dependence on 
angle of attack is illustrated in Figure 15 for the 
smooth airfoil and Figure 16 for the sparsely rough- 
ened airfoil. Figure 15 shows IRT with spray air 
Frossling number data for the smooth airfoil at a 
1.2x10* nominal Reynolds number, for 0°,2°,4° and 
6° angle of attack. Data for the suction side of the 
airfoil are represented by positive s/c values. The 
stagnation gage heat transfer generally increases 
with angle of attack while the first pressure gage 
decreases going from 0° to 2° to 4°, then increases 
from 4° to 6°. This behavior can be explained by the 
movement of the aerodynamic stagnation point; the 
stagnation region sees an effectively larger leading 
edge radius as it moves toward gage 3 with 2° and 
4° angle of attack, and this results in a lower heat 
transfer coefficient. The flow is then accelerated 
around the leading edge, thus increasing the heat 
transfer at the geometric stagnation point. 
Similarly, going to 6° moves the aerodynamic stag- 
nation point toward gage 2, causing an increase in 
gage 3 heat transfer. It should also be noted that 
the power law curve fits showed that as the 
aerodynamic stagnation point moves closer to the 
pressure side gages (gages 2 and 3), the Re°^ depen- 
dence is seen to exist at those gages, i.e., the 
Nusselt number at the aerodynamic stagnation 
point always has a Re 05 dependence. On the suc- 
tion side the heat transfer generally decreases with 
angle. However at 6°, the data illustrate a drastic 
rise in heat transfer at s/c greater than 0.06, which 
is most likely due to boundary layer transition. 

The sparsely roughened airfoil data, mean- 
while, exhibited similar behavior concerning the 
aerodynamic stagnation point, as evidenced by 
Figure 16. The downstream gages, however, 
showed the heat transfer increasing with increasing 
angle-of-attack. The other roughness cases exhibit- 
ed similar angle of attack behavior. 

Data were also acquired at a -4° angle of at- 
tack which allowed more complete measurement of 
pressure side heat transfer values. Comparing the 
-4° (pressure side) data with the +4° (suction side) 
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data showed that for the smooth airfoil the pressure 
side heat transfer was slightly higher than the suc- 
tion side, while for the roughened airfoil the pres- 
sure side experienced lower heat transfer than the 
suction side. 

Experimental Data Versus Cylinder and 

Jlat Plate Correlations 

In computer codes, heat transfer from airfoils 
is often estimated by using cylinder-in -cross-flow 
heat transfer values in the leading edge region and 
flat-plate heat transfer values further aft. Figures 
17 and 18 show chord-based Frossling number 
values corresponding to Frossling s analytical cylin- 
der-in-crossflow solution [16], together with laminar 
and turbulent flat-plate values [17] on the same 
graphs as the present 0°, smooth-airfoil flight and 
tunnel data (without spray air), respectively. Both 
graphs show good agreement with the laminar flat 
plate correlation at s/c greater than 0.06, the flight 
data averaging around 8% lower and the tunnel 
data around 3% higher. However, in the stagnation 
region the experimental heat transfer is somewhat 
lower than that predicted by Frossling. Gages 2-5 
show a 19-21% lower Frossling number for the 
flight data and 14-17% lower for the tunnel data. 
In addition, both the flight and the IRT data fall 
between the laminar and turbulent flat-plate values 
for 0.01 < s/c < 0.06. It seems, therefore, that the 
inscribed cylinder method for estimating the heat 
transfer from the forward portion of a NACA 0012 
airfoil does not accurately work; in this case the 
method substantially overpredicts the measured 
heat transfer. 

Figure 19 shows the IRT data for the dense 2, 
0°, without-spray condition compared with the 
cylinder and flat plate heat transfer values. The 
heat transfer in the stagnation region for the dense 
2 roughened airfoil agrees fairly well with 
Frossling’s smooth-cylinder laminar flow solution. 
Moving downstream on the airfoil, the heat transfer 
drastically increases, reaching a maximum level 
near s/c of 0.035, and then decreases to a level fairly 
consistent with turbulent flow flat-plate heat 
transfer values. The measured heat transfer at 
specific Reynolds numbers are somewhat higher 
than their respective flat- plate turbulent values. 
However the higher measured heat transfer may be 
due to the increase in surface area due to the 
roughness elements (3-7% increase on each gage for 
the dense roughness patterns) that was not taken 
into account in the data analysis. It may be 


mentioned here that the maximum heat transfer is 
in the same general region, if slightly aft, of ice horn 
growth observed during glaze ice accretion [2]. 

Conclusions 

Local heat transfer measurements from a 
smooth and roughened NACA 0012 airfoil were 
successfully obtained in flight and in the NASA 
Lewis Icing Research Tunnel using the method and 
apparatus described in this work. 

M^jor conclusions resulting from this study 

are: 

1. The smooth -airfoil Frossling number data for 
flight with a measured turbulence intensity of 
<0.1% and for the IRT with a 0.5-0.7% turbulence 
intensity showed fairly good agreement at the lower 
Reynolds number (Re=1.2xlO G ). At the higher 
Reynolds number (Re=2.4xl0 6 ), the IRT data were 
somewhat higher than the flight data. Comparison 
of the flight and tunnel rough-surface data showed 
that the effect of turbulence was a slight increase in 
heat transfer. 

2. The IRT is a fairly clean wind tunnel, with a 
measured turbulence intensity around 0.5-0.7%. 
The addition of spray air to the tunnel flow did not 
change the heat transfer, indicating that the 
turbulence level was not significantly altered by the 
spray air. 

3. The addition of roughness to the airfoil 
surface drastically increased the heat transfer 
downstream of stagnation. The roughness elements 
disturbed the laminar boundary layer flow and in 
some cases caused a transition to turbulent flow, 

4. Generally, the suction-side heat transfer on 
the smooth airfoil slowly decreased with angle of at- 
tack, however, at 6° the downstream heat transfer 
drastically increased indicating that some flow 
transition had occurred. The rough ened-airfoil 
cases showed the suction side heat transfer mono- 
tonically increasing with angle of attack. The 
Nusselt number for the smooth-and roughened-air- 
foil cases at the aerodynamic stagnation point al- 
ways correlated with Re 0 5 . 

5. The flight and tunnel smooth-airfoil data 
show good agreement with the laminar flat-plate 
heat transfer values for s/c equal or greater than 
0.06. In the leading edge region, the measured heat 
transfer is somewhat lower than that predicted by 
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Frossling's laminar flow cylinder solution: 15% 
lower in the IRT and 20% lower in flight. 
Therefore, it would appear that the method of using 
an inscribed cylinder for approximating the leading 
edge heat transfer does not work for the NACA 
0012 airfoil. 
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black AND white photograph 



Figure 1. NACA 0012 test airfoil* 


HEAT RUX GAGES -7 



Figure 2. Cross-section of NACA 0012 airfoil with heat flux gages, 


Table 1* Location and surface area of heat transfer gages. 


GAGE# 

* S/G 

SURFACE 



AREA 

sq.cm. 

1 

-0.036 

3.145 

2 

-0.024 

3.145 

3 

-0.012 

3.145 

4 

0 

3.187 

5 

0.012 

3.145 

6 

0.024 

3.145 

7 

0.036 

3.145 

8 

0.048 

3.145 

9 

0.060 

3.145 

10 

0.072 

3,145 

11 

0.083 

3.145 

12 

0.095 

3.145 


2 
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Dense 1 roughness pattern 


Dense 2 roughness pattern 



Typical roughness element 


Figure 3. Location of roughness elements relative to heat flux gages. 
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Figure 4. NASA Lewis Icing Research Tunnel. 



Figure 5. Test airfoil positioned in the Icing Research Tunnel. 
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Frossling Number 



Gage Number 

Figure 9. Frossling Number versus dimensionless surface 

distance; sparse roughness, 0° angle of attack, Re*l-2E6, 
flight and IRT data. 



Gage Number 


Figure 10. Frossling Number versus dimensionless surface 

distance: sparse roughness, 0° angle of attack, Re*2«4E6, 
flight and IRT data. 



Gage Number 


Figure 1 1. Frossling Number versus dimensionless surface 
distance: dense 1 roughness, 0° angle of attack, 
Re-1-2E6, flight and IRT data. 



Gage Number 

Figure 12. Frossling Number versus dimensionless surface 
distance: dense 1 roughness, 0° angle of attack. 
Re «2^4E6, flight and IKT data. 
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Frossling Number Frossling Number 



Gage Number 

Figure 13. Frossling Number versus dimensionless surface 
distance: dense 2 roughness, 0° angle of attack. 
Re »L2E6, IKT data. 



Dimensionless Surface Distance 

23456789 10 11 

Gage Number 

Figure 14. Frossling Number versus dimensionless surface 
distance: dense 2 roughness, 0° angle of attack. 
Re -2AE6, HIT data. 



Gage Number 

Figure 15. Frossling Number versus dimensionless surface 
distance: smooth airfoil. Re- 1.2E6, with spray air, 
IKT data. 



Gage Number 

Figure 16. Frossling Number votsus dimensionless surface 
distance: sparse roughness. Re *24E6, without 
spray air, IRT data. 
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Frossling Number 



Figure 17. Frossling Number comparison for smooth airfoil, (f angle 
of attack: flight data versus Dimensionless Correlations 
for Cylinder [16] and Flat Plate [17]. 



Figure 18. Frossling Number comparison for smooth airfoil, 0° angle 
of attack: IRT data (without spray) versus Dimensionless 
Correlations for Cylinder [16] and Flat Plate [17]. 



-0.04 -0 ,02 0.00 0 02 0.04 0.06 0.08 0.10 

* Dimansionless Surface Distance 

Figure 19. Frossling Number comparison for dense 2 roughness, 9° angle 
of attack: IRT data (without spray) versus Dimensionless 
Correlations for Cylinder [16] and Fiat Plate [17]. 
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